Coulomb drag between two spin incoherent Luttinger liquids 
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In a one dimensional electron gas at low enough density, the magnetic (spin) exchange energy 
J between neighboring electrons is exponentially suppressed relative to the characteristic charge 
energy, the Fermi energy Ef- At non-zero temperature T, the energy hierarchy J <^T <^ Ef can 
be reached, and we refer to this as the spin incoherent Lutinger liquid state. We discuss the Coulomb 
drag between two parallel quantum wires in the spin incoherent regime, as well as the crossover to 
this state from the low temperature regime by using a model of a fluctuating Wigner solid. As the 
temperature increases from zero to above J for a fixed electron density, the 2kF oscillations in the 
density-density correlations are lost. As a result, the temperature dependence of the Coulomb drag 
is dramatically altered and non-monotonic dependence may result. Drag between wires of equal 
and unequal density are discussed, as well as the effects of weak disorder in the wires. We speculate 
that weak disorder may play an important role in extracting information about quantum wires in 
real drag experiments. 

PACS numbers: 71.10.Pm,71.27.+a,73.21.-b 
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I. INTRODUCTION 

In recent years correlated electron systems at the 
nanoscale and in reduced dimensions have attracted 



much attention 



1.2.3 



In one spatial dimension electron 



correlations are expected to be enhanced, leading to the 
so called Luttinger liquid (LL) staten^ The existence of 
the LL liquid state in a one dimensional electron system 
is now an established experimental fact^Si^ with direct 
measurements of the distinct spin and charge velocities 
in momentum resolved tunneling (as predicted by the 
theory) providing compelling evidence 

Another way to explore the correlations in one di- 
mensional systems is in a drag experiment between two 
parallel quantum wires or nanotubesiiiii^ii^ (Recall that 
in most cases the bare conductance of a quantum wire 
is 2e^/h per transverse channel regardless of the elec- 
tron interactions, ^"^'^^ and so does not reveal information 
about electron correlations. An exception to this quanti- 
zation condition is the situation discussed by Matveev in 
Refs. ll^llTi ) The typical drag set-up involves a current 
driven in a "active" wire while a voltage drop is measured 
in a "passive" wire. See Fig. 

The quantity often taken to describe the drag effect is 
the "drag resistivity" (drag per unit length) defined as 



,. 1 dV2 
lim 

h^o h L all 



(1) 



where V2 is a voltage induced in wire 2 (the "passive" 
wire) due to a current Ii in wire 1 (the "active" wire). 
Here e is the charge of the electron, h Planck's constant, 
and L the length of the wire. The sign of the drag can 
be cither positive or negative, but it is generally positive 
(note minus sign in formula) for repulsive interactions 
between electrons. 

Typically, one is interested in the dependence oi 
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FIG. 1: Coulomb drag schematic. Two quantum wires are 
arranged parallel to one another. A dc current 7i flows in 
the "active" wire 1 and a voltage bias V2 is measured in the 
"passive" wire 2 when I2 = 0. 



on the temperature, the interwire spacing, the electron 
density and Fermi wave vector in each wire, the disor- 
der, the wire length, and possibly on an external mag- 
netic field. Physically drag is the result of collisions (mo- 
mentum transfer) from electrons in the active wire which 
tend to "push" or "pull" the electrons in the passive wire. 
Electrons in the passive wire move under these collisions 
until an electric field is built up in the passive wire (due 
to a non-uniform density of electrons there) which just 
cancels the force of the momentum transfer of the elec- 
trons in the active wire. This is the physics behind the 
well known drag formulas of Zheng and MacDonald,— 
and Pustilnik et al^ (See Eq. 0.) 

While the experimental data on drag between quan- 
tum wires is limited jii^iSii^ a fair amount of theoret- 
ical work has been published. Various studies have 
made use of LL theory^SSiSLSSiSiSi Fermi liquid (FL) the- 
ory with^^ and withoulii^ multiple sub-bands, effects of 
inter-wire tunneling)2& effects of disorder^Si shot noise 
correlationSj^^'^™ mesoscopic fiuctuationsj^S*^ and the ef- 
fects of different signs of electron exchange interactions 
in the wires >^ Additionally, phonon mediated drag has 
been studied in one dimension. '^i^ The main qualita- 
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tive difference that is found between the LL and FL ap- 
proaches is whether tends to increase (LL) or decrease 
(FL) as the temperature is reduced to the lowest values. 
For the case of two infinitely long indentical clean wires 
the following results are obtained: In a LL, electrons tend 
to "lock" into a commensurate state at the lowest tem- 
pertures giving rise to a diverging drag, while in a FL the 
phase space available for scattering tends to zero as the 
temperature is lowered, implying a vanishing drag. For 
two non-identical wires the drag is usually significantly 
suppressed at low temperatures relative to the drag in 
the identical caseii2iS£ 

In spite of the theoretical effort, a number of open ques- 
tions remain. In particular, the drag effect is known to be 
strongest when the electron density is lowfii*i2ii^ which 
typically implies that electron interactions are strong, or 



equivalently that = a/(2aB) is large with a 



and 



gb the Bohr radius. For very strong interactions, there is 
an exponential separation of the spin exchange energy, J, 
and the characteristic charge energy, Ep, which at finite 
temperatures can lead to incoherent (thermally excited) 
spin degrees of freedom while the charge degrees of free- 
dom remain approximately coherent and close to their 
ground stateiiSiii This energy hierarchy at finite tem- 
perature T, J<CT<Ci^FjWe refer to as the spin inco- 
herent Luttinger liquid regime. Already there is mount- 
ing understanding of how such spin incoherent Luttinger 
liquids behave in the Green's function^SiSLSI in the mo- 
mentum distribution functiouf^ in momentum resolved 
tunneling, ■^^ and in transport liSiiLiMi 

Our goal in this work is to explore some of the im- 
plications of spin incoherence on the drag between two 
quantum wires. We consider only the simplest case of a 
single channel wire. We attempt to elucidate what qual- 
itative and quantitative changes one can expect for the 
Coulomb drag when the temperature is much smaller or 
larger than J. Since J/Ep is expontially smalljiSiii a 
small change in the temperature can induce a dramatic 
change in the temperature dependence of the drag. Based 
on earlier work of the authors,^*^ we are able to discuss 
the drag deep in the spin incoherent regime in terms of 
spinless electrons using a simple mapping between the 
charge variables of the charge sector of a LL with spin 
and the variables of a spinless LL. 

The crossover to the spin incoherent regime is dis- 
cussed using a model of a fluctuating Wigner solid with 
an antiferromagnetic Heisenberg spin chain in the spin 
sector. Distortions of the solid couple the spin and charge 
degrees of freedom. The model allows us to quantita- 
tively address the crossover. The main result is that 
as the temperature increases from zero to above J for 
a fixed electron density, the (already weak) 2kF oscilla- 
tions in the density-density correlations are rapidly lost. 
As a result, the drag is dramatically suppressed (when 
kpd ^ 1) since forward scattering contributions vanish 
in the LL model and the Wigner solid model, and ikp 
contributions are suppressed relative to 2k p contribu- 
tions by ij{AkF)/ij{2kF) ^ e-'^'""^ for kpd > 1. Here 




FIG. 2: Schematic of the possible temperature dependence 
of the Coulomb drag in a wire of sufficiently low electron 
density that J <^ Ep ^ ujq- The non-monotonic tempera- 
ture dependence shown can be obtained for Kc > 1/2 and 
U{4:kF) <^ U{2kF) which may be realized for widely space 
wires. The temperature T* is the "locking" temperature of 
two identical wires, below which the drag exhibits activated 
behavior and Es is an energy gap associated with the "lock- 
ing". In this paper we consider only temperatures x > T* 
and investigate the limits T < J and T > J, as well as the 
crossover between the two. When J <^ Ef a sharp drop in 
the drag resistance should be observable for T ^ J. 



U{kd) is the Fourier transform of the interwire electron- 
electron interaction, d is the interwire separation, and 
kp = 7r/(2a) is the Fermi wavevector. See Fig.|21 

In addition to a possible dramatic suppression (de- 
pending on various physical parameters) of the drag, a 
non-monotonic dependence on temperature may also oc- 
cur, as illustrated schematically in Fig. |21 Specifically, 
when T > T*, where T* is the "locking" temperature 
(see Sec. im below), we find that the temperature depen- 
dence of the drag is given by 



rD = CiT''''^fiT/J) + C2T° 



(2) 



where f{X ^ 0) ^ 1 and f{X > 1) - e-"^ , with c a 
constant of order one. The coefficients Ci oc U{2kF) and 
C2 oc U{AkF). The exponents a^kp/Akp depend on the 
interactions in the wires and on the presence or absence 
of disorder. We find 



and 



2Kc — 1 clean, 
2Kr disordered 



8Kc — 3 clean, 
8Kc — 2 disordered 



(3) 



(4) 



While it is interesting that disorder changes the tem- 
perature dependence of the drag, it may play an even 
more important role in the measurement of the drag ef- 
fect itself in actual experiments. The reason is that the 
drag effect is maximal for identical wires; any kp mis- 
match results in a drag that is generally strongly dimin- 
ished relative to this case, indeed, exponentially so at 
low temperature. Disorder tends to "smear" the momen- 
tum structure of the density response that determines the 
drag, eliminating this exponential suppression. Hence, 



3 



since in a real experiment one will never have truly iden- 
tical wires, some residual disorder may actually play a 
key role in experimental studies of drag between one di- 
mensional systems. 

This paper is organized in the following way. In SecHTI 
we discuss general considerations for drag between two 
quantum wires with an emphasis on features relevant to 
the effects of being in the spin incoherent regime. In 
Sec, mil we discuss a model of a fluctuating Wigner sohd 
with a Heisenberg spin chain to describe the magnetic 
(spin) degrees of freedom. We derive expressions for 
the density fluctuations of the electron gas by includ- 
ing magneto-elastic coupling that induces 2kF density 
modulations at low temperatures (T <C J) and results 
in an instability towards a local lattice distortion favor- 
ing a spin-Peierls-like state. In Sec. IIVI we discuss the 
temperature dependence of the Fourier transform of the 
density-density correlation function in detail. In Sec. 
we discuss the drag itself in detail by considering different 
temperature regimes, the effects of disorder, and the case 
of wires with mismatched density. In Sec. lVII we summa- 
rize the main results of our paper and in the appendicies 
we give some exact formulas for the density-density cor- 
relation function relevant to the spin coherent-incoherent 
crossover as well as other useful formulas. 



II. GENERAL CONSIDERATIONS 

We assume the Hamiltonian of our system is of the 
form 



H = Hi+H2 + H: 



12 



(5) 



where Hi is the Hamiltonian in the i*^ wire and H12 
describes the interactions between electrons in different 
wires. A proper drag situation is one in which the tun- 
neling between wires can be neglected. We thus assume 
that H12 allows only interactions which forbid electrons 
to tunnel between the wires4S The Hamiltonian Hi in 
principle describes arbitrary interactions between elec- 
trons within the i*'' wire; depending on the particu- 
lar situation of interest, a number of different models 
have been proposed from Fermi liquid^iiSiS^ to Luttinger 
Iiquid»-Sm4 

The Hamiltonian H12 is a function of the interwire 
electron interaction which is often taken to be of the form 



Ul2{x) 



(6) 



where e is the charge of the electron, e the dielec- 
tric constant of the material, and d the separation be- 
tween the two wires. The Fourier transform, U{k) = 
2{e^ / e)KQ{kd), depends on the dimensionless parameter 
kd and for kd ^ 1, lJ{k) ~ -^=e~'''^. This exponen- 
tial dependence of the interwire interaction leads to an 
exponential dependence of the drag resistivity on wire 



separation (for large enough d) when the drag is dom- 
inated by large momentum transfer as it is in the LL 
model and the model we will discuss in this paper. The 
following drag formula (or its equivalent) has been de- 
rived by Zheng and MacDonald^^ for a disordered FLj^ 
by Klesse and SternSi for a LL, and by Pustilnik et alm^ 
for a clean FL: 



dk 



"° ^^ k^UUk) Imxf(fc,c.)Imx?(fc,c.) 



sinh^(tj/2r) 



(7) 



where Imxf^(A;,w) is the imaginary part of the Fourier 
transformed retarded density-density correlation func- 
tion, Eqs. (|33|l and (|34|l . and rii is the density of the elec- 
trons in the i*'' wire. Knowing the general features of the 
interwire interaction t/12, (see Fig. it is clear that to 
determine the drag one must determine the Imxf (fc,w). 
We will discuss the behavior of Imxf (fc, w) and its tem- 




FIG. 3: Momentum dependence of the the quantity 

k^Ui2ik) = (^^) ixKo{x))'^ appearing in the drag formula 

Eq. Q as a function of the dimensionless variable x — kd 
where we have assumed the real space inter-wire electron 
interaction @. Plotted is the quantity {xKo{x))^ which 
has a maximum for x ~ 0.6. For x <^ 1, {xKo{x))'^ ~ 
(— 7-|-ln(2) — ln(a;))'^x^ where 7 f» 0.5772 is Euler's constant, 
while for a; > 1, {xKo(x)f ~ e~^^. 

perature dependence in detail in Sec. IIVI but for now it 
is sufficient to point out that most of weight occurs near 
momenta of A; « 0, fc « 2fci?, and k « -ikp so that one 
may write 



Imxf (fc, uj) « Imxf '°(A:, uj) + Imx"'^"^ {k, uj) 

+lmxf'''^'"'{k,u;). 



R,2kf 



(8) 



As we will see, the 2kF components of Imxf (A:, ^) present 
at low temperatures, T <^ J, will disappear when T 3> 
J. Moreover, once it is known that lmxf''''\k,uj) cx 
S{u; — Vc,ik) (see Sec. IIV|) where Vc.i is the charge ve- 
locity of the J*'' wire, it is readily seen from Eq. ^ that 
for the LL model (or any model possessing a harmonic 
bosonized theory) the k = contribution to the drag 
resistivity vanishes (in the absence of disorder), leaving 
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only contributions from the higher k = 2kF,^kF val- 
ues of Imxf (fc, w). The higher k — 2kp,^kp values can 
be strongly suppressed by the interwire interaction be- 
cause C/i2(fc) ~ when fee? ^ 1. Hence, in the limit 
kpd ^ 1 we expect a dramatic reduction in the drag 
and a change in the temperature dependence of the drag 
(see below) when 2fci? oscillations are lost and only Akp 
modulations remain. However, before we address these 
details, it is useful to go over what we can say about drag 
deep in the spin incoherent regime itself. 

Based on the results of earlier workf22i2i we can imme- 
diatly discuss the drag deep within the spin incoherent 
LL regime, J <C T <C Ep- It has earher been arguedift 
that the spin incoherent LL behaves essentially as a spin- 
less LL by noting that the Hamiltonian is diagonal in spin 
to lowest order in J /T <^ 1 , and by demanding the equiv- 
alence of the physical charge density in both cases. The 
equivalence can be summarized by the simple equation 
K = 2Kc (in the notation of Klesse and Stern^^) relating 
the interaction parameter of the spinless LL theory and 
the interaction parameter of the charge sector of the LL 
theory with spin. (Strictly speaking, for drag the relevant 
interaction parameter is K = K- , the interaction param- 
eter in the odd channel of the relative density of the two 
wireS)22i2i but when the inter-wire interactions are suffi- 
ciently weak K- is essentially equal to the correponding 
parameter of the isolated wire.) The relation K — 2Kc 
is valid for any particle conserving operator;i2i the drag 
resistivity is derived from such an operator. Hence, those 
general results apply her e. 

As discussed in Refs. I20l2ll drag in the spinless LL 
model comes from the backscattering of electrons. As 
long as the inter-wire interactions are sufficiently small 
the effects of backscattering can be treated perturba- 
tively. Based on such a perturbative treatment, Klesse 
and Stern foundSi that for identical wires there is a tem- 
perature scale T* that separates high and low tempera- 
ture drag regimes. (It is assumed througout this paper 
that T* is greater than the thermal length Tl so that 
the Fermi liquid leads attached to the quantum wire are 
not felt.) For temperatures T ^ T*,Tl, the drag varies 
as a power of the temperature, while for temperatures 
T* ^ T ^ Tl, the drag resistivity shows activated be- 
havior with a gap of order T* itself. The physics is simi- 
lar to that of a pinned charge density wave. This result 
follows from an analysis of a sine-Gordon model in the 
odd channel of the coupled wire problem. Applying the 
equivalence rule K = 2Kc discussed above in the spin 
incoherent regime, we have 

p^^r8^=-3^ t:$>t*,Tl,j, (9) 



PD-e^^^^, t*->t:^TlJ, (10) 

where E'^ ~ T*,^ Note that for 3/8 < < 3/4 the 
temperature dependence of the spinless (fully polarized) 
electron gas and the spin incoherent electron gas exhibit 



very different drag resistivity behavior with temperature 
when T » T* . (The spinless case has a diverging drag 
resistivity as T is lowered, while the spin incoherent case 
has a suppressed drag resistivity as T is lowered.) This 
is qualitatively similar to the transport results found in 
Ref.EOfor 1/2 <gc<l. 

The results (O and H1U|) above were derived from a 
perturbative analysis of the sine-Gordon equation which 
results from treating the backscattering in LL theory. 
For most realistic parameter values, the baskcattering 
strength flows to strong coupling and the resulting state 
is that of the two quantum wires locked into a "zig- 
zag" charge pattern. The value of T* depends on de- 
tails of the quantum wire system such as the density, 
wire widths, and separation cZjSi but for most realistic 
situations T*/Ep < 0.01. 

It is interesting to consider how spin incoherence af- 
fects the "zig-zag" locking pattern of the electrons in the 
two wires. The relative size of J and T* will determine 
what the periodicity of the "zig-zag" pattern will be for 
T « T*. For J < T < T*, there is a Hkp" lock- 
ing (seen easily from the K — 2Kc mapping^) since 
T ^ J ensures 2k p pieces of the density are washed 
out, while for T < J,T*, there is a "2^^" locking. Of 
course, for T* ^ J <C T the locking phase is not ob- 
tained. Throughout this paper we will assume T ^ T* 
so that we need not be concerned with "locking" from 
here forward. 

Similar arguments to those given in Ref. |23 can also 
be used to describe the incommensurate-commensurate 
transition deep in the spin incoherent regime for wires 
of different electron densities. We now leave generalities 
behind and turn to a detailed calculation of the drag itself 
in the regime of very strongly interacting one dimensional 
electrons. 



III. THE FLUCTUATING WIGNER SOLID 
MODEL 

We assume from the outset that the interactions be- 
tween the electrons are very strong, which typically 
means the density is low enough that we can treat the 
electrons in each wire as a harmonic chainM in the charge 
sector and a nearest neighbor Heisenberg antiferromag- 
net in the spin sectortii^ 

i?wire — He + Hs , (11) 

where 

^2 2 

^^^^l^^^rr^ (12) 
1=1 

is the Hamiltonian in the charge sector with pi the mo- 
mentum of the electron, ui the displacement from 
equilibrium of the electron, m the electron mass, and 
luq the frequency of local electron displacements (this will 
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depend on the electron density, the width of the wires, 
the dielectric constant of the material, and other param- 
eters such as the distance to a nearby gate^^*^). The 
position of the electrons along the chain are given by 



xi = la + ui , 



(13) 



where a is the mean spacing of the electrons. The Hamil- 
tonian of the spin sector takes the form 



i+i 



(14) 



Note that in H14|l the coupling Ji between spins depends 
on the distance between them. Assuming that the fluctu- 
ations from the equilibrium positions are small compared 
to the mean particle spacing, we can expand the exchange 
energy as 



Ji ^ Jo + Ji(m/+i - "0 + 0{{ui+i - uiY) . 
In this case the full Hamiltonian takes the form 



where 



H, 



+1 



ui)Si 



i+i 



Si 



(15) 



(16) 



(17) 



Here H^-c represents a magneto-elastic coupling as it 
couples the magnetic modes to the elastic distortions of 
the lattice that constitute the charge modes. 

Our goal is to evaluate the Fourier transform of the 
retarded density-density correlation function —i9{t — 
t'){[p{x,t), p{x' ,t')]) [which appears in the drag formula 
O] up to second order in Ji for T, Jq <C hcoo, for both 
T <C Jo and T ^ Jq. We use the following definition 
of the electron density: p{x,t) = J2iH^ — al — ui{t)). 
An exact calculation within this model is presented in 
Appendix ^ Here we will pursue an approximate cal- 
culation that captures all of the essential features of the 
more exact perturbative results. 



A. Low energy approach to charge fluctuations 

In this work we are concerned only with energies (tem- 
peratures) small compared to the characteristic lattice 
energy, i.e. T <C Tiujq, but still large compared to the 
"locking" temperature T* . When T ^ Jq further ap- 
proximations that can be made, but for now our only 
restriction will be that T <ti htOQ. We begin by expand- 
ing the displacement of the electron density in a Fourier 
series. For low energy distortions the k component 
is most important, while the magneto-elastic term H17|l 
couples the k n/a component to the spin operator 
Si+i ■ Si. Thus, the displacement, ui of the Z*'* electron 
in the harmonic chain l|12|l is approximately given by 



ui = uo{la) + U7r(Za)(— 1)' 



(18) 



where uq refers the fc « component of the displacement 
and refers to the k ^ ir/a displacement. Both uq and 
Utt are assumed to be slowly varying functions of position, 
and we expect u-^ <C uo- 



1. Low energy form of the action 




FIG. 4: Approximate form for the phonon spectrum. There 
are two relevant parts of the spectrum that enter the action 
Sc, the part near fc ~ and the part near k ~ n/a. The 
phonons with fc « are described by a Debye model while 
those phonons with k ~ n/a are described by an Einstein 
model. 



The action for the low density electron gas is 

s = Sc Ss -\- Ss—c ■ 



(19) 



Using the expression H18|) for the particle displacement 
fields, and noting that from the phonon dispersion of 
(|12|l . ujk = 2a;o| sin(A;a/2)|, one has for fc « the dis- 
persion iVk — Vc\k\ with Vc = ujQa^ while for k « n/a the 
dispersion cuk — 2ujo is independent of k. (See Fig. 0]) 
In the charge sector we then have the following action of 
the form Sc = S^ 



Q7T 



Sc^ 



dxdr 

2lT 



KcVc 



{{d,ec{x,T)f + vl{dJc{x,T)f 



{{drUTr{x, r))^ + {2u!o)'^u^{x, rf 



(20) 



where Kc = ttTi/ {2amvc). Note the lack of spatial deriva- 
tive in the u^r piece of the action which results in the 
absence of any k dependence in ujk near k n /a. In ef- 
fect, we have described the small k oscillations (phonons) 
with a Debye-type model and the large k oscillations as 
an Einstein model. Fig. ^illustrates the approximations 
to the full phonon spectrum. 

In the standard LL model for weakly interacting elec- 
trons Kc = vf/vc ^ 1 where is the Fermi velocity of 
non-interacting electrons. In the present case of strongly 
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interacting electrons we have Kr = „^ % y-!— ~ ~ 

~ 7- where E'^ is the Fermi energy of non-interacting 

electrons and = 0/(205) where ub = eh^/me'^ is the 
Bohr radius of a material of dielectric constant e. Thus, 
in the strongly interacting limit Kc scales roughly as the 
ratio of the kinetic energy to the potential energy which 
itself roughly scales as rj^ implying that very strong 
(long-range) interactions can lead to small Kc- In prac- 
tice, however, Kc rarely appears to be smaller than 0.2 
or so. 

By making the identification £1 = (— l)'S'/+i • Si 
£{x,t) in the continuum limit, we have spin-charge cou- 
pling in the action 



Ss — c 



dxdr 
27r 



2JiUt^{x, t)£(x, t) 



(21) 



In the special situation where T <^ Jq, the action for 
the spin sector (|14(l can be bosonized a o^Ti"^? 



dxdr 

271 



K.v, 



{idr0s{x,T))' +v!{d,e,ix,T)f 



(22) 

where the SU{2) symmetry of the Heisenberg model 
implies that Kg — 1 and the spin velocity is Vs ~ 
Joa/h. (Here intrawire backscattering effects in the 
spin sector have been neglected and a sine-Gordon term 
dropped. Since we are ultimately interested in ener- 
gies/temperatures much larger than Jg or much smaller 
than Jo this will not affect any of our conclusions.) How- 
ever, when T ^ Jq the action for the spin sector must 
describe more accurately the short distance physics of 
the Heisenberg chain. Nevertheless, the action H22|l will 
prove useful in understanding the approach to the spin 
incoherent regime in the limit T — > Jq from below. 

Finally, for T ^ Jq the coupling term in the action 
can be expressed as^ 



Ss-c^ I ^^7^u^(x,r)sin(V20,(x,r)) , (23) 
Ztt Znaa 



where we have used the low-energy bosonized form 
i-iySi+i ■ Si « ^sm{V2es{x)). Here a = 0{a) is 
a short distance cut off of order the lattice spacing. At 
low energies this term will lead to the spin-Peierls state 
indicated in Fig. |5| 



2. Expressing the density 

Expanding the density and making use of (|18|l then 
j;ives (we have suppressed the time dependence of uq and 



UjT immediately below for clarity of presentation) 



p{x) = 



I 



I 



I 



S{x ~ la — ui), 

S{x — la ~ UQ{la) + U7r(^o)(— 1)'), 
[S{x — la — uo{la)) 

~S'{x -la- uo{la))u^{la){-iy)]. (24) 



Multiplying by d{x' — la), integrating over x', and using 
the Poisson summation identity. 



S{x' -la) = - 



(25) 



m— — 00 



the expression for the density becomes 
p(x) « / dx' —5{x — x' + uo{x')) 

m—~: 

X 1 + COs{ — x')dx'UTr{x') 

{{2m + l)e*^^' -f (2m - l)e 
2o V 



a 



u^{x') , (26) 



where we have integrated by parts in the second term of 
(|24|l . Performing the integration over x', making use of 
the approximate relation 8(x — x' — uo{x')) « S{x' — x — 
uo{x))/\l + dxUo{x)\, and assuming adxUTr{x) <C Ut^{x), 
we find the most important terms up to Akp are 



p{x) w -(1 - dxUo{x)) 
a 



2tt 



Po 



— sin (2fci?(a; + uq{x))) Ut^Ix) 
a 

+ cos{AkF{x + Uq{x))) , 

—dx9c{x) — pq- — sin (2kpx + \/29c{x) \ Ut;{x) 
TT o V / 

+/O0 cos {^^kpx 4- VS9c{x)^ , (27) 



V2 



where po = kp = n/(2a), and we have made the 
identification UQ{x)/a = \/2Qc[x)I'k. Recall the field 9c 
is governed by the action 12U|I . This formula resembles 
the standard bosonized expression for the density of a 
Luttinger liquid (except for the term u-^ multiplying the 
2kF part of the density instead of a term involving the 
spin fields). As we will now see, a formula very simi- 
lar to that obtained from the standard Luttinger liquid 
treatment results from integrating out the high energy 
phonon modes in favor of the low energy spin modes. 
To see this, consider only the 2kp part of the density and 
compute {p2kp{x,T)) to lowest order in the action Ss~c 
and integrate out the fields to obtain a new p'^^ (cc, t) 
independent of u^r. At lowest order we find 



(u^(x,t)) 



2Ji 



dx'dr' 
2tt 



£{x',r') 



x{TrUT,{x,T)u^{x\T'))i 



(28) 
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where Tr is the r-ordering operator and the r-ordered 
product is evaluated in the action given in (|^. The 
r-ordered product is readily evaluated as 



(u^(x,T)u^(a;',r'))s; 

g2uo|T-r'| 



hS{x 



X 



P-2ujq\t-t'\ 



as /3a;o 



(29) 



Recall that we are interested only in temperatures low 
compared to the phonon energy wq so the limit /Swo oo 
is the appropriate one. Here /? = {kBT)~^ where ks is 
Boltzmann's constant. The integral over position in H28|) 
is immediately evaluated with the delta function in H29|l 
and the remaining integral over r' can be approximately 
evaluated under the assumption Ploq oo which pro- 
duces the dominant contribution at r' « r with a width 
in r' of order 1/{2llIq) resulting in 



2Ji 



£ix,T) 



(30) 



which yields 



(x,r) 



1 



sm 



{2kFX + V2ec[x,T)^ 

x8{x,t). (31) 



The result (|31|) is general, and valid whenever T <C 
htoo. However, when T Jq one may use the expression 
£{x, t) = sin(\/26's(a;, r)) which leads to the familiar 
looking density 



P {x,t) = po 



—dx9c{x,T) 

TT 



IGtt 



Ji 



sin ( 2kFX + \/26'c(2;, r)j sin(\/26's(x, r)) 

-l-po cos (^4fcFa; + ^/S9c{x, t)^ . (32) 



that when the temperature is low compared to both spin 
and charge energies a 1-d electron gas always behaves 
as a LL in the sense of the various power laws that will 
appear in the correlation functions, although the overall 
prefactors of the 2kp pieces will be down by the ratio 

(rruf^) ■ other hand, the system is in the 

spin incoherent regime Jq <^ T ^ Ticjo, the 2k p parts of 
the correlations will be washed out from thermal effects. 
We now turn to an investigation of how this happens in 
detail for the case of the density-density correlation func- 
tion and then discuss the implications for the Coulomb 
drag between two quantum wires. 



IV. EVALUATION OF lmxfik,uj) 

Here we consider two limits of the double wire system 
shown in Fig. ^ (i) Clean wires without disorder and 
(ii) Wires with weak disorder. The case of strong disor- 
der is uninteresting as the electrons are all localized over 
the relevant energy/length scales of the experiment. As 
we discussed in Sec. |h1 the drag formula Q generically 
contains contributions at fc « 0, fc ~ 2kp, and k « Akp 
so that lmxf{k,uj) w Imxf"'^ {k , lu) + Imxf^''^'^'' {k,uj) + 
Imxi ' {k,uj). We now turn to an evaluation of each of 
these pieces. We have used the two standard (equivalent) 
definitions 

/oo /"OO 
dx / die*«'^+"')*-'=^) 
-00 Jo 

X ([(p,(x, t) - p,,o), {p^{0, 0) - p,,o)]), (33) 

and 

/oo 
dx / dre^('^"^-'=^) 
00 Jo 

X {{p,{x, t) - p^,o)(MO, 0) - p,,o)), (34) 



The expression for the effective density H32f) with the 
high energy modes integrated out in favor of the spin 
variables is valid for T ^ Jq and may be compared di- 
rectly with the LL result obtained for weakly interacting 
electronsi^ At temperatures T 3> Jo (|ST|l must be used for 
the 2kp density variations. The only material difference 
between (|32f) and the standard LL result is the dimen- 

sionless ratio of spin and charge energies (^ ^^i^2 ^ ~ ^ 
which is absent (since it is of order 1) in the familiar LL 
case. When the interactions are strong as we have as- 
sumed them to be here, then ( — ) ~ — <C 1, since 

Ji diminishes and uiq increases with increasing strength 
of the interactions. However, starting from the strongly 
interacting limit and decreasing the interaction strength 

-> 1. It is worth emphasizing, then. 



the ratio 



Ji 
— 



where pi,o is the average density of the i'^ wire, and rj 
is a small infinitesimal that ensures convergence of the 
time integral in (|33|l . The retarded correlation function 
is obtained from (|34|) via the substitution za;„ —^uJ + ir|. 
We will use both of the formulas above in the subsections 
that follow. 



A. Clean Wires 

We first consider wires with no disorder. We will 
also assume initially that T ^ Jq so that we may use 
the form of the density (|32|l . (This is only an issue 
for the evaluation of Imx^'^*^^ (fc, w) since Imx f' ''^ {k,uj) 
and Imxf ''^'^^ {k, oj) do not involve the spin sector of the 
Hamiltonian.) As the authors discussed in Ref. ^3 the 
approach to the spin incoherent regime from tempera- 
tures well below the spin energy can be understood in 
this way. In all calculations below, recall that we have 
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assumed the temperature is low, T <C Tiwo, so that the 
charge sector is always in the LL regime and described 
by the action .9° in if^ . 



Imxf'°(fc,a.) 



We first evaluate the fc « piece of the retarded 
density-density correlation function. From the expres- 
sion we have 

pf{x,t)^pa-—d^e,{x,t), (35) 

TT 

whose correlation function is readily computed (see Ap- 
pendix 0) to yield 



Imxf'*'(fc,w) 



fc2 fiL 
a? 2muJk 



[■k5{lj — bJk) — nS{uj + LUk)] ■ 



(36) 

The equation above, H3()l) . is the central result of this sub- 
section and it is worth pausing to emphasize some of its 
features. Most notably, while the calculation was done at 
finite temperature, there is no temperature dependence 
of Imxf'''^ (kjOj). Thus, the finite temperature fc « re- 
sponse is identical to the zero temperature response. This 
means that temperature does not "broaden" the zero 
temperature 5-function reponse. Moreover, for the model 
at hand, at small |fc|, ujk = Wi.cl^l so that for a given k 
there is a unique value of Uk- This means then when the 
result (|36|) is substituted into the drag formula {Tj) the 
drag is identically zero. (An exception is the measure 
zero point where the wires are identical, i.e. vi^c = V2,c, 
and the drag response is infinite. For real wires this pre- 
cise matching is not possible and the k = part of the 
drag generically vanishes.) 

Ultimately, the vanishing of the fc w drag is a result 
of the delta functions appearing in It is expected 

that the delta functions will be broadcnediSi^ in a more 
complete treatment and that this will lead to a non-zero 
and temperature dependent fc w contribution to the 
drag. 

In our work here, we have assumed from the outset 
that the electron interactions are very strong and a di- 
rect bosonization of the electron operator is not valid. 
Instead, the approximation we have made to obtain the 
action 12(J|) . which is formally identical to that obtained 
for weakly interacting electrons with a linear disper- 
sion (aside from the Ujr terms), is to treat the displace- 
ments of electrons to lowest order in the Taylor series: 
{ui+i — ui)/a ~ \/2dx0c{x)/'n. Including higher deriva- 
tives would result in an interacting bosonic theory and 
would likewise broaden the delta functions in (|36|l by 
an amount inversely proportional to the lifetime and 
would yield a finite fc « drag. The precise nature of 
this contribution to the drag is still a subject of ongo- 
ing research>i2i^ It is therefore difficult to compare it 
quantitatively in theoretical calculations to the 2kp and 
Akp contributions. However, we expect that it may be 



larger or smaller than the latter depending upon circum- 
stances. For instance, the fc « drag is clearly sub- 
dominant for drag between identical, clean wires, with 
repulsive interactions. Fortunately, for our purposes of 
discerning the spin coherent to incoherent crossover at 
T w J, we may satisfy ourselves with the observation 
that the fc w drag is in any case featureless at this tem- 
perature. Hence, it can easily be "subtracted" by looking 
for strong temperature-dependent changes in the drag in 
this temperature window. 



The 2kp component of the density response and its 
temperature dependence is the central issue in this pa- 
per and we now turn to it in detail. We have already 
discussed general features of the spin incoherent limit 
T Jo in Sec. |nl and we will discuss other more de- 
tailed and quantitative features of that regime in the next 
subsection where we consider Imx^'^''^ (fc, w). Here, we 
will initially assume that the temperature is low com- 
pared to the spin energies, T <^ Jq, and use the low 
energy density expression H32|l . Starting from the low 
temperature limit we show that as the temperature be- 
comes of order the spin energy, the temperature depen- 
dence of the 2kp part of the drag changes and rapidly 
vanishes as Jq ^ for fixed T ^ Jq. We also show 
that in the low temperature limit we recover the tem- 
perature dependence of the drag obtained by Klesse and 
Stern^i for electrons with spin. When kpd ^ 1 the loss 
of 2k p contributions to the drag (when T Jq) implies 
(via Eq. j?)) and Fig. O that there is expected to be a 
dramatic reduction in the drag over a very small temper- 
ature window when only the Akp contribution remains, 
as U{Akp)/ij{2kp) - 6-2'=^'' for kpd > 1. 

The 2k p part of the low energy density operator H32|) 



Ji 



Po 



sin {2kpx + V20c{x, r)^ 

X sin(%/26l,(a;,r)) (37) 



which leads to the following finite temperature result for 
the 2k p part of the density-density correlation function 
computed from and 



-z([p^t(x,t),p^X(0,0)]) = (^ 



Ji 



X cos(2fc^2;)Im 



(TrTa/Wf 



.[sinh (^(x - «ei)) sinh (^(a; + v.t)) 

(Trra/t;,)^- 

[sinh (^{x - Vst)) sinh (^{x + Vst))]^ 



Here a is a short-distance cut off of order the lattice 
spacing. We note that in Eq. (|38|) - and in subsequent 



similar formulae - singularities at x = ±Vct, ivst are reg- 
ularized by infinitesimal imaginary parts to the time t, 
which for ease of presentation are not shown. It is worth 
pointing out that because of the hyperbolic nature of the 
correlation function at finite temperature, a temperature 
dependent "coherence length" naturally appears in both 
the spin and charge sectors. From inspection, the charge 
coherence length ^dT) = Vc/{KcTrT) ~ ahujo/{kBT), 
and the spin coherence length ^s(T) = Vs/{KsnT) ^ 
aJ^/ (ksT). Strong interactions imply Vg/vc ^ 1 (Jo 
hujo) so that ^^(T) < £,c{T). Note that ^^(T) « a when 
r« Jo. 

Our task is now to substitute (|38|l into the integral in 
and evaluate the integrals over position and time. 
Unfortunately, this integral does not appear to have a 
closed, analytical form. Nevertheless, its general struc- 
ture is apparent. At zero temperature the structure in 
the {k,Lu) plane is very similar to that of the Green's 
function already computed by Voit^a and by Meden and 
Schonhammeri^ Depending on the value of there are 
singularities or thresholds at w = Vsk± and u = Vck±, 
where k± = k ± 2kp. With small, but finite temper- 
ature these features are smoothed out. However, as 
the temperature increases towards Jo the overall weight 
'^'^'^^(fc±, begins to rapidly diminish. To see 

Jo from below. Then 



this, consider the limit T 
xf'^'°^(fc±,a;) can be bounded as 



xf— )" 
V167ry 



keT 
Jo 
Ji 



oo 

dx 



xlm 



sinh - M)) sinh (^ix + M))]"^ 



(39) 



where 

R,2kF 



2kf 



cx 



Jl 



{2kFfUl^{2kF)T 



(40) 

where we have used the result that at low enough 
temperatures, dkk'^Ul2{k)lmxi{k,ijj)lmx2{k,ijj) ~ 



\ ?f ^ ^ ^ 



c is a constant of order unity. For fixed T, 
Xi '~ ' {k±,^) ^ as Jo —> 0. This conclusion is in- 
dependent of the particular form of the operator used 
in the spin sector. For example, using the more general 
expression H31f) will lead to the same conclusion for any 

£{x, t). Thus, the already weak (because {^ ^'^1^2 ^ ^1) xlm 
2kp density oscillations are rapidly suppressed with tem- 
peratures once T ^ Jo- See Fig. |31 for an illustration. 

Having emphasized how "fragile" xf'^^'' {k± , w) is for 
T ~ Jo, let us now return to the low temperature 
limit T ^ Jq. In this limit, the temperature depen- 
dence of Imxf'^''''(fc,w) = x^''^''''{k+,uj)+xf''^''''{k-,i^) 
can be readily extracted by making the substitutions 
X — ttTx/vc and t = nTt and then computing the Fourier 
transform. With this substitution, we find 



G) G> G> 



FIG. 5: Classical Wigner solid model. The electrons are 
indicated by solid black dots and are separated by a mean 
spacing distance a. Top: A "cartoon" of the antiferromag- 
netic spin arrangement is indicated by the arrows. The lattice 
has density oscillations of smallest wavevector ikp- Bottom: 
Magneto-elastic coupling allows the system to lower its energy 
by slightly distoring the lattice in order to gain magnetic en- 
ergy. Stronger spin correlations are indicated by the ovals 
which pair neighboring electrons at spacing less than a. The 
lattice has density oscillations of smallest wavevector 2fcF, in- 
dicating twice the period of the undistorted lattice. When 
T ^ Jo the (small) lattice distortion is thermally washed out 
leaving only the 4A;_f periodicity of the underlying Wigner 
solid. 



{2kF)^Ui2{2kF)linxi^{2kF,uj)lmxf{2kF,uj)Ak, where 
Afc ~ T and used the result that the oj integration in 
(O for l|38() does not contribute to any temperature 
dependence of the drag. The function f{X^O) 1 
and f{X > 1) ~ e^'^^ , where c is a constant of order 
unity. 

The result (|40|l is identical to the result obtained by 
Stern and KlesseSi in the weakly interacting limit of the 
1-d electron gas when T <C Jo- Note that while the 
temperature dependence is the same in the low temper- 
ature limit, the overall result is still down by a factor 

~ I — I ^ 1 when the interactions are strong. 

For completeness, it is worth emphasizing that in the 
high temperature regime (T ^ Jq) the expression H31|l 
must be used for the 2kF part of the density. In this case, 
one must compute the Fourier transform of the correlator 



^{[pTkA^,t),pf,^io,o)]) 



Jl 



16 V moj^a^ 



cos{2kFx) 



{TiTa/vc) 



sinh (^(x - v,t)^ sinh (^(a; -t- v,t) 



x([£:(x,i),f(0,0)]) 



(41) 



In the high temperature regime (|41|l will not behave much 
differently from H38|) when T w Jo. In particular, we 
expect 



:[f(:E,0),£:(0,0)]) 



(42) 



where ^ a and so in the high temperature limit 
the results will be qualitatively similar to what we dis- 
cussed earlier. Of course, the detailed structure of 
Imxf '^'^^ (fc, w) for T w Jo requires that (|41|l be used. 
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This in turn requires that the dimer-dimer correlation 
function {[£{x,t),£{0,0)]) be evaluated by a more gen- 
eral (perhaps numerical) method than the effective low 
energy theory given in (|22|l . 



3. Imxf'*'=-(fc,^) 



to obtain 



RAkp/i \ 



p|a^ / 27rTa y^-"' r(l - 



47rT 



-)nKc 



■ uj — Vck± 

^ 4:7TT ' 



r(i-if,-zi£+£^)r(i 



(48) 



4ttT 



In the previous subsection we saw that when v^/vc <^ 
1 and temperature T ^ Jq, the 2kp contributions to 
the drag are dramatically suppressed and only the -ikp 
contributions remain. In contrast to the case of the 2kp 
density fluctuations, the present model H20|l allows for a 

closed analytic expression for Im^^'^*^^ (fc, ti^)- We begin 
with the 4kp part of the density operator 132|) , 

(4kFX + V8e,ix,t)y (43) 

which leads, after evaluating the correlators at finite tem- 
perature, to 



xlm 



i{[plt t), Pit (0, 0)]) = pI cosiAkpx) 



sinh 



(^(x - Wei)) sinh (^(x + i;,t)) 



2Ka 



(44) 



As in the case of Imxf^''^'"' {k,uj) the temperature de- 
pendence at low enough temperatures can be extracted 
by making the substitutions x = ttTx/vc and t = irTt. 

{k,uj) 



This then leads us to x 



RAkp/i \ 



Xt (k±,uj)=T 



x4lm 



dx / dte 

oc Jo 



[sinh(x — i) sinh(a; + t)] 



, (45) 



and k± = k ± Akp. By the same arguments made in the 
previous subsection (that the form substituted into 
Q leads to no temperature dependence of the drag from 
the uj integration, and that the dominant contribution 
from the k integral comes from k± with Afc ^ T) 
the temperature dependence of the Akp contribution to 
the drag is 



r%^^ (X {4kprU't^i'ikp)r 



8K^ 



(46) 



which is identical to the result Q obtained in Sec. HTIbv 
applying the general arguments of Ref. 40 for the map- 
ping of a spin incoherent LL to a spinless LL. 

Fortunately, the Fourier transform (|45|l can be com- 
puted exactlyiSiS This is done by making the change of 
variables si = x — t and S2 — x + t, and using the integral 
result^ 



ds- 



sinh(s)]s 



25- 



_, T{g/2-iz/2)T{l-g) 

r(i-.9/2) 



(47) 
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FIG. 6: (Color online) Frequency dependence of 

Imxf ''"'^ (fc± = l,!^), Eq. 14811 . for various interaction 
strengths and temperatures T. The charge velocity is 
fixed at He = 1. At Kc = 0.5 there is a crossover from a sharp 
peak for Kc < 0.5 a.t u = k± to monotonically increasing 
(with Lu) threshold-type behavior for Kc > 0.5. Finite tem- 
perature acts to smear the T — lj = k± singularity and adds 
weight to Imxf'^''^(fc± = l,u>) for lo < k±. 

The temperature dependence of Imx^'^'^'^ (A:± = 
for different interaction values Kc is shown in Fig. O At 
Kc = 0.5 there is a crossover from a divergence to a 
threshold-type behavior. The main effect of the temper- 
ature is to smooth the sharper features near co — Vck. 



B. Weakly Disordered Wires 

1. Slowly varying background potential 

Modulation doping in quantum wire systems gives rise 
to a smoothly varing background potential. Such disor- 
der has an important effect on the drag as it impacts the 
nature of the electronic states that participate in dragi^ 
Since the coupling of the density to the potential de- 
pends crucially on the Fourier components k, there is 
an important difference in how the charge density cou- 
ples to disorder in the spin coherent and spin incoherent 
regimes. Consider the following coupling of the density 
to background potential modulations: 



Ha 



dx 



dxV{x)p{x) 

Vo{x)po{x) + V2feF (x)p2feF (a;) + V'4A;f (a;)p4feF (2^) 



(49) 
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Here V2kF {x) 



iAkpx'\ 



i2kFX'\ 



and Vikp {x) 



We can study the scaling dimensions of 
V2fej7, Va^Up using the expression for the density, Eq. 
after integrating out the high energy field u-j^ in favor 
of the lower energy spin fields. The scaling dimensions 
of the different scattering terms can then be determined 
from the action (where the integration over x has already 
been caried out) 



'dis 



dTV2 



2kj 



and 



dTV. 



which gives 



dim[V2kp] = 1 - 



h.c. 



El 

2 



dim[V4fc^] = l-2Kc 



(50) 



(51) 



(52) 



(53) 



In these units SU{2) invariance implies Kg = I, so that 
the 2k F piece is more relevant than the 4:kF piece of the 
potential whenever Kc > 1/3. Thus, we expect to see 
strong temperature dependence of the pinning of the den- 
sity whenever Kc > 1/3 as the more relevant 2fci? piece 
wiU be lost for T > Jq. Moreover, if 1/2 < < 1 the 
2kp piece is relevant while the Akp piece is irrelevant. In 
this case, the effect should be most dramatic. The regime 
1/2 < Kc < 1 can be reached for large but finite [/ in a 
one dimensional Hubbard model. 



2. Effects of random correlations in forward scattering on 
the correlation functions 

As the fc w part of the back ground potential fluctu- 
ations are often the most important at low energies, it is 
worthile to review^ their influence on the density-density 
correlation function. The forward scattering part of the 
background potential is 



V2 

dxVo{x)po{x) = / dxVo{x)d^6c{x) , 



(54) 

where we have we have used the result ()32f) and assumed 
J dxVo{x)po = 0. Forward scattering can be completely 
eliminated from the Hamiltonian (action) by making the 
change of variables 



0cix) = edx) 



V2Kc 



dzVo{z) 



(55) 



and completing the square in Eg. 1)20(1 . Assuming that 
the disorder has white noise correlations given by the 
distribution Pv„ = exp[— I?^^ J dz|Vo(z)P], 



D , 



Vo{x)Vo{x') ^ -5{x - x') 



(56) 



where the overbar indicates a disorder average, and the 
constant D Uc/7'scatt with Tscatt the typical scattering 
time for the electrons. The disorder averaged parts of 
the density-density correlation function can then readily 
be determined: 



{[pf{x,tlpf(QM) = {[pf{x,t),pf{QM)\vo=o 



{[pfkpM^pfkpim])- 

e-^(^)'l-l([p§t(x,t),p°t(0,0)])|y„.o, 



(57) 



(58) 



{[pt, (x, t),pll (0, 0)]) I vo=o. (59) 

It is evident that larger wavevectors are suppressed more 
by the forward scattering with no suppression at all for 
the k Ri part of the density. Treating the 2k p and 
4:kp backscattering contributions with white noise corre- 
lations analogous to (|56|l is more involved and requires 
studying renormalization group flows. However, we re- 
iterate that if the disorder is slowly-varying (as expected 
from donor potential modulations in modulation doped 
structures), the 2kp and Akp potentials are relatively 
weak and probably negligible, at least in the sorts of 
structures optimized for the drag measurements we envi- 
sion here! 



3. Fourier transform of disorder averaged correlation 
functions 

Once the Fourier transforms of the 2k p and ikp 
density-density correlation functions l|38|) and (|44|l are 
known, the Fourier transforms of the disorder averaged 
correlation functions are readily computed from the con- 
volution theorem 



Imxf (fc,t^) 



dq 



21 



_^ 2^l-^ + {k-qY 



Imxf(<z,c^), (60) 



where l^^ — D{Kclvc)^ for the 2kp pieces and l^^ — 
AD{Kc/vc)'^ for the ^kp pieces. The main effect of the 
disorder is thus to broaden the singularities in x,f^(fc,w) 
by an amount of order in momentum space. Hence, 
the 4fci? singularities are broadened 4 times as much as 
the 2kp singularities. 

Fig. [3 illustrates the effect of weak disorder on 
Imxf '^'^^^ (fc, (jj). Qualitatively the effects of disorder are 
very similar to finite temperature-there is a broadening 
of the sharpest features in the response. This has impli- 
cations for the temperature dependence of the drag as we 
discuss in the next section. 



V. STUDY OF THE COULOMB DRAG 

We have already discussed several features of the drag 
in the previous sections of this paper, including the tem- 
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FIG. 7: Disorder and thermal broadening of Imxf '^''^ {k, ui = 
1). The charge velocity is fixed at = 1. The disorder broad- 
ening is computed using Eq. . The efltects of weak forward 
scattering on the 4A;f density fluctuations are qualitatively 
similar to that of thermal broadening-see case of / = 0.06 kp- 
Analagous results are also obtained for the 2k f density fluc- 
tuations. 



perature dependence (in certain limits) of the k 2kp 
and k « Akp contributions to the drag. Implicit in 
those discussions was that the wires were identical. In 
this section we will present numerical calculations of the 
Coulomb drag as a function of temperature for identi- 
cal and non-identical wires, and attempt to illuminate 
the crossover to the spin incoherent regime with semi- 
quantitative estimates. 



Drag at low temperatures and in the crossover 
to the spin incoherent regime 

1. The low temperature Luttinger liquid regime 



From here on, we will assume SU{2) symmetry which im- 
plies Kg = 1- When Kg — 1 it is clear that the 2kF con- 
tribution to the drag is an increasing function of temper- 
ature whenever Kc > 1/2 and a decreasing function oth- 
erwise. For the Akp contribution the boundary between 
increasing and decreasing contributions is Kc — 3/8. Fi- 
nally, by comparing the exponents of the 2kp and Akp 
terms, one finds that the Akp pieces dominate the drag 
for Kc > 1/3 when T > T** , while the 2k p pieces dom- 
inate the drag for Kc < 1/3 in the same temperature 
regime. Note that this implies that the 2kF density os- 
cillations are more important for the drag at higher tem- 
peratures when the interactions are strong enough that 
Kc < 1/3. This requires, of course, that the system is 
still at low enough temperatures that the spin degrees 
of freedom can be described by the effective low energy 
theory H22() . In order to obtain T** < Jq and for the 
analysis above to be reasonable, we expect that we must 
have C/(2fcF) > ij{4.kp). 

From the results of Appendix El we can express the 
ratio 



B{Kc,Ks,Vs/vc) 
A{Kc,Vc) 



Ji 



l2kF{Kc,Ks,Vs/Vc) 



■ {Kc 



(63) 

where l2kF{Kc, Ks,Vs/vc) and hkpiKc) are given by 
Eqs. (|C1|I and HC2I) . As Vg/vc the crossover 
temperature (|62(l becomes very small because both 

l2kAKc,Ks,Vs/vc) < hkAKc) and ( ig^^^Lga^ ) ^ 

(t)^) ^ ^ thai limit. Of course, as Vg shrinks for 
fixed Vc, the temperature range over which the LL the- 
ory itself is valid is also shrinking and the spin incoherent 
regime (where only the Akp density modulations remain) 
is approached. 



At the lowest temperatures where T <^ Jq we showed 
that the low energy theory H2U|) . H22|l results in the fol- 
lowing temperature dependence of the drag resistivity Q 



TD ^ A{Kc,Vc){^kpfuU^kF) 
+B{Kc,Kg,Vs/vc){2kFfuW2kp) 



( T 



2r~r2 --• - ' ^ 



8/fe-3 



(61) 



where A{Kc, Vc) and B(Kc, Ks, Vg/vc) are functions that 
depend on the variables indicated and T^j^^ = ira / ^JVcVg 
and T^/ — ira/vc are effective temperatures in the re- 
spective sectors. It is evident from (161(1 that the temper- 
ature dependence of the k w 2kp and k w 4fcF contribu- 
tions are different so there is a temperature at which the 
two balance out: 



2. Crossover to the spin incoherent regime 

The hallmark of the spin incoherent regime is the 
equivalence of the real electron system with spin to 
a spinless system.^" with the exception of the Green's 
function^ and other non-particle conserving operators. 
In the case of drag, spin incoherence manifests itself as 
a thermal washing out of the 2kp oscillations in the 
density-density correlation function (|SJ|. When the in- 
teractions are as strong as they are here, the weight 
of the 2k p oscillations are already down by a factor 

2 

even at zero temperature. 



T* 



T2kf 



4kf 



\ T2kf 



8Ka 



B{Kc,Kg,Vg/vc)Uh{2kp) 
AAiKcVcM^i'^kp 



(62) 



As we have discussed before )2i the spin incoherent 
regime can be understood by starting with T 3> Jo, tak- 
ing Jo ^ 0, for fixed T and then finally taking T 0. 
In the present formulation this is equivalent to fixing a 
jfinite, but low temperature, applying the low energy the- 
ory (|20|) and H22|) and then taking the limit Vg/vc — > 
as we did in the previous section. The approach to the 
spin incoherent drag regime can be directly obtained via 
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this procedure. One expects that as Vg is lowered, the 
power law H4UI) will first breakdown (at temperatures it 
once held for larger v^) before the contribution vanishes 
altogether from /(T/Jq). 



B. Drag in the spin incoherent regime 

In this subsection we present some numerical results 
justifying earlier analytical arguments for the tempera- 
ture dependence of the drag. We first consider identical 
wires and then we study non-identical wires. 



O 



O 



. _ = 0.3. Eq.(4) 

— K =0.3.rpaT'"'-"' 

■ ■ ■ K. = 0.3. Eq.(4), 1"' =0.13kp 



- K =0.6, Eq.(4) 
- K =0.6,rj, a t"*- 
■■ K = 0.6, Eq.(4),f' =0.13 kp 




-3 -2 

Log(T) 



1. Identical Wires 

When the wires are identical we expect the tempera- 
ture dependence of the Coulomb drag given in Eq. l|HT|l 
to be obtained at the lowest energies. However, as we 
have seen in the previous sections the 2kp oscillations 
are rapidly washed out in the limit Vg/vc —> and only 
the Akp oscillations remain. In this subsection, we pro- 
vide numerical evidence that the manipulations leading 
to the Akp temperature dependence of (|61|l is justified. 
Since these are also the same arguments leading to the 
2kp temperature dependence of r/^ at the lowest tem- 
perature, these are implicitly justified as well. Fig. |S1 
illustrates the comparison between the exact result from 
(|^ substituted into and the approximate power 
law (|46|l . A disorder value of = 0.13^^ was used 
in Eq. (|60|l to compute the drag of the disordered system 
from Eq. lO . The drag was computed over a temperature 
range O.OIwq < T < O.Swq- For fc^T > hvc/l the drag 
of the clean and disordered system are indistinguishable, 
while for fc^T <C TiVc/l there is a crossover of the temper- 
ature dependence to another power of the temperature. 
Empirically, we found the power law 



4fc_F, disorder rp%Ka — 2 



(64) 



to be a very good fit for any value of < ifc < 1- This 
temperature dependence can actually be inferred from 
(|45|l . lO and Fig. EJ As we have argued several times 
earlier, the lo integration in |(7J) does not contribute any 
temperature dependence beyond the (7^4/s:c-2-j2 factors 
in front of H45I) (with the square coming from the drag 
formula When /csT <C hvc/l, the k integration 

picks up a contribution proportional to rather than T. 
Adding up the exponents leads to 



.disorder ^ rp8Kc 



FIG. 8: Temperature dependence of the Coulomb drag in 
the spin incoherent regime where only the 4kF components of 
the density fluctuations contribute. The drag was computed 
over a temperature range O.Olojo < T < O.Sojq. Shown is a 
comparison between the drag formula Q with H48|l substi- 
tuted in and the approximate power law 14611 . Two values of 
the interaction parameter Kc are shown. A disorder value of 
l^^ = O.lSkp was used in Eq. 1601 to compute the drag from 
Eq. lO. For fcsT ^ hvc/l the drag for the clean and dis- 
ordered system are indistinguishable, while for ksT <C hvc/l 
there is a crossover of the temperature dependence to another 
power of the temperature. Empirically, we found the power 

to be a very good fit for any value 



law rtf ^ 



disorder ^ rp8I<c 



oiO<Kc<l. 



the spin incoherent regime. In Fig.|51we present some nu- 
merical results for the dependence of the drag for small 
density mismatches between the two wires. Note that 
weaker interactions and higher temperatures lead to a 
more robust drag effect between two wires of slightly dif- 
ferent densities. Note also that with only a few percent 
change in the relative densities of the wires the drag effect 
is substantially reduced. 

The effects of disorder on the drag for density mis- 
matched wires is shown in Fig. ^| When Hvc/l ^ ksT 
the disorder has a significant effect-making the drag more 
robust for non-identical wires. While for small \ni — 712! 
the drag is reduced relative to the clean limit, for larger 
\ni ~ 712! there can be appreciable enhancement. 

Finally, we note that in the case of clean wires of dif- 
ferent charge velocities but the same Fermi wavevector 
(in which case the Kc are diffferent) the temperature de- 
pendence of the drag is 



F ^ y4(K,,i+if,,2)-3 



2. Drag for non-identical wires 



and if the temperature is also much less than Jg, 



' D 



(65) 



(66) 



Coulomb drag for non-identical wires and the 
incommensurate-commensurate transition has been dis- 
cussed for fully coherent clean wires in Ref. Via 
the mapping detailed in Ref. ^3 the incommensurate- 
commensurate transition can be ready discussed deep in 



VI. SUMMARY OF RESULTS 

We have discussed the Coulomb drag between two 
quantum wires in the limit of low electron density where 
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FIG. 9: (Color online) Coulomb drag dependence on relative 
wire density at fixed temperature. We have set Vc = 1. The 
Coulomb drag between two wires drops rapidly as a function 
of density mismatch. Only a few percent change in the rel- 
ative densities of the wires results in a dramatic suppression 
of the drag. At higher temperatures and weaker interactions 
the drag response is more robust to small density differences 
between the wires. 
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FIG. 10; Coulomb drag dependence on relative wire density 
at fixed temperature for different values of disorder. Here 
T = 0.02 Wo and Vc = 1. While for clean wires the Coulomb 
drag drops rapidly as a function of density mismatch, some 
amount of residual disorder allows for a more robust drag 
effect between quantum wires that may have slightly different 
values of fcF- 

at finite temperatures the energy hierarchy J <^ T <^ Ep 
can be obtainediS In this limit, the spin degrees of free- 
dom are completely incoherent and we have shown this 
implies the loss oi2kp oscillations in the density-density 
correlation function. As a result, a non-monotonic tem- 
perature dependence of the drag on temperature may re- 
sult. In the spin incoherent Luttinger liquid regime, the 
drag problem maps onto the identical problem for a spin- 
less Luttinger liquid only with Kc^ — > 'i.Kc- so that for 
a clean wire the drag resistivity goes as po oc 7"8i<'c--3^ 
where Kc- is the coupling parameter of the antisymmet- 
ric charge mode of a Luttinger liquid theory with spin. 
We have shown this temperature dependence explicitly 
with approximate analytical calculations and confirmed 
those approximations with numerics. 



Our results are based on a fluctuating Wigner solid 
model appropriate to quantum wires in a very strongly 
interacting regime, which typically implies low electron 
density. The spin sector is modeled as a nearest-neighbor 
antiferromagnetic Heisenberg spin chain. Without any 
coupling of the spin and charge sectors, no 2k p den- 
sity modulations appear. However, including a magneto- 
elastic coupling term that allows for a linear change in 
the nearest-neighbor spin coupling for small distortions 
induces 2kp oscillations in the density. The coupling is 
weak, however, and the 2kp oscillations are easily washed 
out by temperatures T ^ J . 

The fluctuating Wigner solid model is studied by de- 
riving effective expressions for the density operator when 
the highest energy phonon modes are integrated out in 
favor of spin operators. At the lowest energies (includ- 
ing the spin energy) an expression is obtained equivalent 
to that known well in Luttinger liquid theory except the 
2kp terms contain a prefactor of order the ratio of ki- 
netic to potential energy. Nevertheless, in this limit all 
correlations functions exhibit power law decay with the 
familiar exponents of the spin and charge sectors. 

These density operators are then used to compute 
density-density correlation functions which are Fourier 
transformed into frequency and momentum space and 
used in previously derived drag formulas. Since the den- 
sity operator has contributions at momenta fc « 0, few 
2k and k Akp the Coulomb drag will generically con- 
tain contributions from each of these pieces. We show ex- 
plicitly that the A; « piece vanishes due to the harmonic 
approximation to the Wigner solid. This is equivalent to 
linearizing the electron dispersion in the standard Lut- 
tinger liquid treatment for weakly interacting electrons. 
Generically, the k 2kp and k Akp contributions are 
non-vanishing and we explicitly compute their tempera- 
ture dependence, t'^^c-Iq-cT/Jq g^^^ ySKc-a^ ^j^g ^ow 

temperature regimes. 

We have also considered the case of non-ideal wires in 
which weak disorder is present. We find that white-noise 
correlated forward scattering disorder does not affect the 
k result, while it tends to broaden the sharp fc— space 
features of the k « 2kp and k k, Akp density-density cor- 
relation function in a manner similar to temperature. As 
a result, the drag resistivity crosses over to a different 
power law, T^Xc g^j^j which is increase by one 

power of the temperature relative to the clean result. 
Finally, we have also studied the reduction of the drag 
due to a density mismatch between the two wires and 
shown the drag may be substantially reduced with only 
a few percent change in the relative densities of the wires. 
When disorder is present the drag is more robust to den- 
sity mismatches between the two wires and this fact is 
likely to play an important role in real drag experiments. 

We hope that this work will help to inspire further 
experimental studies on one-dimensional drag, which to 
date is quite limited. 
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APPENDIX A: EXACT EXPRESSIONS FOR 

{p{x,r)p{x',r')) UP TO SECOND ORDER IN Ji 

The low-energy description given in Sec. IIIII can be 
treated more accurately, but in a less physically trans- 
parent way by applying the results of this appendix. 



1. Diagonalization of He 

The charge Hamiltonian H12I) is diagonalized by the 
transformation 



(Al) 
(A2) 



(assuming periodic boundary conditions) where the the- 
ory is quantized by imposing = ihSim and 
[ufciPfc'] — iWkk' ■ Making these substitutions we find 



The general expression for the density-density correla- 
tion function at finite temperature T — 1//3 is 

(p(x,r)p(x',r')) = ^U^P)PMp{^'.r')), ^^^^ 



{umo 



where 



(A9) 

with the operators taken in the interaction representa- 
tion. The averages (...)o = Tr[e~^^°...] where the 0*'' 
order Hamiltonian is (|16|l taken with Ji = 0: 



Hq — Hc + Jo ^ • Si , 



(AlO) 



and H' is the correction to this to be treated in pertur- 
bation theory 



H' = Ji + i - Ui)Si+i ■ Si 



(All) 



3. Evaluation of the density-density correlation 
function 



k—ki 



PkP-k , muj^ 

—Z 1 —UkU^k 

2m 2 



(A3) 



where uo^ — 2aJo| sin (-y-) | ~ Vc\k\ for small k with 
Vc — LOoo. the sound velocity of the charge modes. In 
momentum space the harmonic chain is just a sum of 
harmonic oscillators with frequencies that depend on the 
wavenumber k. 

The Hamiltonian l|A3p can be brought into a particu- 
larly simple form via the transformation 



4. - 

which brings He to 



2n 



niLUk I I 

Uk H P-k 



mujk 



2h 



muJk I 1 

u-k Pk 



rriLOk 



(A4) 
(A5) 



He = ^^fiUk [rik + I 



For later reference it is useful to note that 
ui{t) = Mk{l) (akc-^"^ + al.e^"^ 



(A7) 



where Mk{l) ^ y 



h Akal 



e and we have used ojk 



a. Zeroth Order 

At zeroth order we have Ji = and there is no cou- 
pling between the charge and the spin degrees of freedom. 
Therefore we can completely ingore the spin sector since 
it traces out trivially. Thus, 

{p{x,T)p{x',T'))("^ = Z-^TT[e-^"'p{x,T)p{x',T')] , 

(A12) 

where Zc = Tr[e f^^"] is the partition function of the 
charge sector and the density is expressed as p{x, t) = 
J2i=i 5{x — al ~ ui{t)). From the Hamiltonian (|A6|) Zc 
can be readily evaluated as 



HE 



Za^ w y e 

k nfc=0 



n- 



(A6) which we will make use of later. Thus, 



1 _ e^Phuk ' 
(A13) 



{p{x,T)p{x' y))^°^ = z-^Y. I 

i.i' 



27r 27r 



^iri{x~al) ^i^{x' -al') 



(A14) 



where the quantities ui appearing in the exponent of the 
trace can be expressed using Eq. (jA7|) . Using e^^^ = 



16 



e"^e^e I-^'^l/^ (where [A, B] commutes with both A and 
B separately), we focus on the trace and obtain 

X e-';eS,|M,|^e""('-'')-"5(--') 



g-/3Rwi; (aj.afc + 1/2) g-iC* (/c)o+ g-iC(fe)afc 



(A15) 

where C(fc) = IMkHrje''"'^-'^''''^ ^^ikal'-o,kr')y Evalu- 
ating the trace for each k independently and using defi- 
nition of a Laguerre polynomial of order Uk , 

Lnd-\C{k)f) = (nfc|e-*^*('=)«^e-^^W«'=|"fc) . (A16) 
and then applying the important formula 

1 



cut off when the argument of the In becomes 0(1) which 
occurs when n = ri; ~ (al)"^ + {vc{t — t')Y / {(iUvc)- 
In the limit of T — > 0, /3 — > oo and the second terms 
drops out all together. In this paper we are interested 
in the limit T ^ Ep^ so the second term can be ignored 
altogether. We will not explicitly consider finite temper- 
atures in the first and second order expressions. 



h. First order 

The manipulations needed here are identical to those 
used to compute the zeroth order result, so we simply 
quote the result: 



n=0 

we find 



= -l_exp||Cp^| , (A17) 



ft |M(fc)|^4(l - cos(A;a)) 

27ra 



fe>o 



x(l - cos (6'/c(n,0;r, T')))e 



LI' 



2-K 2-K 



(A23) 



where 



1 + 



1 _ Q-l3h0Jq 



(A18) 



(A19) 



where I is arbitrary, A'' is the number of electrons in the 
system, and 



3=0 



{2j)\{2N-2j)\ \2AJ 



and 



1 _ g,-i3nuiq 



(A20) 

and eq{l, t, f) = qa{l - I') - u!q(t - t'), where t = -it. 
Finally, shifting the summation variables 1 = 1 — 1' and 
performing the integrations we obtain 

{p{x,t)p{x' 



2n]/ F -G{1) 



TT \ (x — x' — al) 

exp 



(27-1)!! 1 

^ A-^-^ , (A24) 

where the n dependence enters through B = i{x — x' — an) 
and A = F — G{n). It is worth noting that neither 
t)p(x',t'))(°' nor {p{x,T)p{x' ,t'))^^^ contain a 2fcir 
component. This component will only appear in the sec- 
ond order term, as we now discuss. 



c. Second order 



The second order corrections are (where x = {x,t)) 



4(F-G(0) i' 
(A21) 



(p(x)p(x'))(^) = Tdn rdT2 
Jo Jo 



where 



F-G{1) = 



ah 



mvr2n 



^a^ + {aiy + ivc{T-T')) 



In 



-(p(x)p(x'))(°) / dn r dT2 

Jo Jo 



[H'(n)ir(r2)p(x)p(ji'))o 

ZcZs 

" {H'{n)H'{T2)}o 



ZcZs 



+2£ln 



n=0 



(g + (n + l)pnva)^ + («0^ + jvcir - r'))^ 
(a + (n + l)mvc)^ 



(A22) 



(r2))(°) 

{H',{l,T^)H',{l\T2)p{^)p{^))0 



Here a is a short distance cut off of order the lattice spac- 
ing a. Note that for finite temperatures the second sum is 



-(P(x)p(x')) 



/N\(0) {H'c{l,T'l)H'^{l' ,T2))o 



(A25) 
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where H'^{1,t) = ui+i{t) — ui{t) is the charge part of 
H' . From (|A25|I it is clear that when dimer-dimer cor- 
relations (S';+i • Si{ti)Si>j^i ■ Si'{t2)Y°'' are present (pre- 
sumably when T ^ Jq) then a 2kp component appears in 
{p{x,t)p(x' ,t')). After some algebra, we reach the final 
form 

1,1' -^0 -^0 



x(Si+i ■ Si{ti)Si>+i ■ Si'{t2)) 



(0) 



/ 27r27r' 



dr] ijj(x-an) i${x' -am) 



xe 

fe 



fc>o 



X Y[ |Mfe|44(l-cos(A:a))2e-2(-i— 2)'^-'^ 

femax 



fc'>0 



(A26) 



where 

h{k) = lAffepe-^^^'''^'' 



R{2n) 



+2riS,\cos[ka{l~l')j 



1 — cos(fca) 



1 



— cos(/sa) 



with 



1 — cos(fca) 
n = n + iojut/ (ka), 
ifi = m -\- ioj^t' / {ka), 



(A27) 

(A28) 
(A29) 



and 

R{s) = cofi{ka{2l' + 2- s)^+ cos{ka{2l' + 1 - s)^ 

-|-cos(^fca(2r-s)y (A30) 

APPENDIX B: COMPUTING lmxf'°(k,Lo) 

Imx^'°(fc, Lo) is computed by making use of the Fourier 
expansion of p'^{x,t), Eq. and the formula 

Consider first the Fourier transform to momentum space: 



dxe-'^- {pf [x, r)pS« (0, 0)) = {pf {k, T)pf (-fc, 0)> , 

where translational invariance was used. The Fourier de- 
composition of p'^{k^T) is readily obtained by making 
use of Eq. H35|) . the relation UQ{x)/a = \/29c{x)/it, and 
the representation of uq{x) given in Ea. IA7I 



where we have implicitly converted the discrete k sums 
to integrals and L is the length of the system. It 
is easily verified that this has the right units to give 
{pf{k,T)pf{~k,Q)) in IHJ the correct dimensions of 
inverse length. Only the expectation values of the cross 
terms a^a^ and al^a_fc are nonzero, giving 



{pf{k,T)pf{-k,0)) 



fc2 hL 



a? 2mL0k 

X [e-'^>'^{l + nB{k))+e^''^nB{-k)\ (B3) 

where nsik) — (afcOfc) — (e'^"'-' — l)^""^ is the boson oc- 
cupation factor. Returning to the expression (|34|l and 
evaluating the t integral we find 



X°(fc,w„) = 



k^ hL 

2mLj}^ 



-1 



1 



(B4) 



which upon the analytic continuation to real frequencies 
iojn Lo + if] leads directly to Eq. (|36|l . 



APPENDIX C: EXPRESSIONS FOR 

l2k^(Kc,Ks,Vs/vc) AND hkp(Kc) 

The functions l2kFiKc,Ks,Vs/vc) and likpiKc) de- 
fined in Eq. I|fi3(l are given by 



hkpiKc^Ks^Vs/vc) 



°° dujE2kAKc,Ks,Vs/vc,^/T) 
T sm\?{oj/2T) 

(CI) 



and 



T sinh2(tj/2r) ' 



(C2) 



where 



OO /'OO 



^2kAKc.Ks,Vs/v,,u/T) = ( Ini / di dte'^' 

-OO Jo 



xlm 



[sinh(5; — i) sinh(a; + t)]' 



[sinh [xvc/vs — t) sinh(5;Wc/ws + i)] 



(C3) 



and 



OO /"OO 



^AkAKc.oj/T) - Im / di die'—' 

\ J-oo Jo 



xlm 



[sinh(a; — t) sinh(a; -I- 1)] 



2K^ 



(C4) 
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